Angle resolved photoemission spectroscopy of Sr 2 Cu02Cl2 — a revisit 
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We have investigated the lowest binding-energy electronic structure of the model cuprate 
Sr2Cu02Cl2 using angle resolved photoemission spectroscopy (ARPES). Our data from about 80 
cleavages of Sr2Cu02Cl2 single crystals give a comprehensive, self-consistent picture of the nature 
of the first electron-removal state in this model undoped CuC>2-plane cuprate. Firstly, we show a 
strong dependence on the polarization of the excitation light which is understandable in the context 
of the matrix element governing the photoemission process, which gives a state with the symme- 
try of a Zhang-Rice singlet. Secondly, the strong, oscillatory dependence of the intensity of the 
Zhang-Rice singlet on the exciting photon-energy is shown to be consistent with interference effects 
connected with the periodicity of the crystal structure in the crystallographic c-direction. Thirdly, 
we measured the dispersion of the first electron-removal states along V — >(tt,tt) and T ^(ir,0), the 
latter being controversial in the literature, and have shown that the data are best fitted using an 
extended t- J-model, and extract the relevant model parameters. An analysis of the spectral weight 
of the first ionization states for different excitation energies within the approach used by Leung et 
al. (Phys. Rev. B 56, 6320 (1997)) results in a strongly photon-energy dependent ratio between the 
coherent and incoherent spectral weight. The possible reasons for this observation and its physical 
implications are discussed. 
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I. INTRODUCTION 

Extensive effort continues to be expended to under- 
stand the electronic properties of the cuprate high T c su- 
perconducting compounds. Apart from many differences, 
the cuprate superconducting materials have important 
properties in common. Namely they are all layered com- 
pounds possessing Cu02-planes, where the origin of the 
superconducting behavior is to be found. The electronic 
structure of the states with the lowest binding energy is 
satisfactorily understandable only in models which take 
into account correlation effects. ARPES has proven to 
be a powerful tool to investigate the low binding-energy 
electronic structure of the cuprates. Important examples 
for the success of this method are the mapping of the 



Fermi surfaced and the determination of an anisotropy of 
the superconducting gap consistent with a d-wave order 
parameterB. 

In the undoped regime the cuprates are two- 
dimensional antiferromagnetic insulators with Cu02- 
planes resulting in a correlation-gapped half-filled 
Cu3dj;2 _ y 2-02p x i '02p y antibonding band of lowest 
binding-energy. Upon doping with holes, the states with 
lowest binding energy evolve towards the chemical poten- 
tial, and the cuprate becomes metallic or superconduct- 
ing. To clarify the main principles of the evolution of 
the lowest binding-energy states with doping by apply- 
ing model Hamiltonians, one has to go step by step from 
simple situations or limiting cases to more complicated 
scenarios. A natural starting point, then, is to investi- 
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gate the propagation of a single hole in a Cu0 2 -plane 
by studying the undoped parent compounds of cuprate 
superconductors by photoemission, where one electron is 
removed from the Cu0 2 -plane. The states thus investi- 
gated can be termed the first electron-removal states. 

A special class of parent compounds are the oxy- 
halides such as Sr 2 Cu0 2 Cl 2 , Ca 2 Cu0 2 Cl2, Sr 2 Cu0 2 F 2 
(one layer compounds), and Ca3Cu 2 04Cl 2 (double lay- 
ered compound)tru, all of which contain an apical halo- 
gen atom, rather than an oxygen. It is now well es- 
tablished that the apical oxygen which is contained in 
most of the cuprate superconductors is not necessary for 
high T c superconductivity as Hiroi et al. ^ showed that 
Ca 2 _a;Na a Cu0 2 Cl 2 is superconducting (T c = 26K). Out 
of the class of cuprate parent compounds without apical 
oxygen, the layer compound Sr 2 Cu0 2 Cl 2 is remarkable 
due to its extremely stable stoichiometry. Despite consid- 
erable efforts so far it has not been possible to dope this 
substance chemically. Furthermore, the Cu0 2 -planes are 
unbuckled and the crystal structure shows neither or- 
thorhombic distortion nor a superstructure at least down 
to 10 K (Ref. g). Therefore, Sr 2 Cu0 2 Cl 2 can be seen as 
the best realization of a two-dimensional antiferromagnet 
at half-filling and it is in this limit an ideal test case for 
any Hamiltonian in the low-doping regime. 

In theoretical studies, the lowest binding-energy ion- 
ization states of the cuprates were predicted to be es- 
sentially described by a singlet antibonding combination 
of a Cu3d 3 .2_j / 2 orbital (containing the intrinsic Cu3d 9 
hole) and a coherent combination of the four neighboring 
02p x /02p y orbitals which can be thought of aSpContain- 
ing the hole created in the photoemission processEm. This 
two-hole state is generally referred to as the Zhang-Rice 
singlet state! . The singlet character of the first electron- 
removal states, at least in CuO, has been experimentally 
verified using spin-resolved resonant photoemissionlij. 

There are various possibilities to theoretically model 
the dynamics of this lowest binding-energy excitation. 
For instance, in the framework of a three-band Hub- 
bard Hamiltonian (H^), the Zhang- Rice singlet state 
is a two-particle state and belongs to the A\ g (totally 
symmetric) irreducible representation of the eigenstates 
of H^b. In the one-band Hubbard Hamiltonian (Hib), 
the t- J-model and its extensions, the Zhang-Rice singlet 
state is an one-particle state belonging to the B\ g irre- 
ducible representation and has the same symmetry as 
a Cx&d x 2 _ y 2 orbital. Figure |l| shows a schematic rep- 
resentation of a Zhang-Rice singlet state, including the 
two mirror planes Mi and M 2 which are perpendicular to 
the Cu0 2 -planc and which are relevant for photoemission 
along along T — >(7r, it) and T 0), respectively. 

It is now generally accepted that the is a good 
starting point to describe the low binding-energy elec- 
tronic structure in the strongly correlated cuprate sys- 
tems and—consequently the has been investigated 
in detailtil. Additionally, it has been shown that the 
i/ib still carries enough information to describe the low- 
energy dynamics in the cuprateala. In the next level of 



simplifying these models, the Hubbard Hamiltonian can 
be expressed in terms of new Fermionic operators. In 
the large-?7-liniit, this new Hamiltonian can be written 
in a t 2 /U series and gives an effective one-band Hamil- 
tonian, the t- J-model. This has been shown to give the 
same results for both the (Ref. 0) and the Hit, (Ref. 
|l3|) approaches. Extensions of the t- J-model have been 
carried out mainly in two different ways. Firstly, two 
further parameters - t' describing the diagonal hopping 
and t" for the next nearest neighbor hopping - can be 
introduced in addition to the direct hopping t. 

The fact that a i-mediated hopping creates strings 
of spin defects leads to a predicted bandwidth for the 
Zhang- Rice singlet state which is governed by t 2 /U ~ J. 
Hopping events mediated by t' and t" , however, take 
place on the same spin sub-lattice and therefore even 
taking small values of these parameters has a signifi- 
cant effect on the results the t- J-model gives for the 
k-dependent spectral function of the Zhang- Rice singlet 
state. Secondly, the t- J-model has also been extended to 
take thcee-site hopping terms (proportional to J/4) into 
accountE-j, which appear in a natural way if one employs 
a more detailed derivation of an effective Hamiltonian for 
the low binding-energy excitationstj. 

A different ansatz to understand the dynamics of holes 
in Cu0 2 -planes was proposed by LaughlinE£l. He argued 
that the hole created by photoemission decays into spin 
and charge degrees of freedom and that the photoemis.- 
sion experiment measures the dispersion of the spinonli3. 
However, in a gauge-field treatment of this model this 
decay of the quasiparticles is suppressed below the Neel- 
temperature, T/v, due to confinement!!-! Another funda- 
mental approach to calculate the properties of the Hub- 
bard Hamiltonian and the t- J-models uses the fact that 
they share an approximate SO (5) symmetryEil 

The preceding discussion highlights the intense the- 
oretical interest in the lowest binding-energy electron- 
removal states of the undoped cuprates. Although orig- 
inally trailing a few years behind the theoretical work, 
experimental investigations of the first electron-removal 
state in|-uridpped cuprates have been carried out using 
ARPESliHci Before going on to mention the status in 
the field to date, we first ask the question: what kind of 
information does an ARPES experiment provide us with? 

Neglecting "extrinsic" effects such as scattering of the 
photoelectrons on their way to the surface and assum- 
ing the applicability of the sudden approximation, the 
photoemission intensity /(k, E), reads 

j(k,E)~ i m *j i 2 s (^ E )f( E ) a) 

\i),\f) 

where the sum runs over all final states and initial 
states, f(E) denotes the Fermi function and S(k, E) the 
spectral function: 

S(k,g) = --Im(»|c k — 1— - 4 |i) . (2) 
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Mi,f — (f | A • p | i) is the matrix element to be 
taken between the initial state | i) with energy Ei and 
final states | /} with energy E, A-p is the photoemission 
interaction operator and H the Hamiltonian of the sys- 
tem. From an ARPES experiment one can gain mainly 
two kinds of information. Firstly, the spectral function 
5(k, E) gives direct information about the dispersion and 
the quasiparticle character of the states and can be com- 
pared directly to predictions of model Hamiltonians. All 
information carried by S(k, E) is expected to be inde- 
pendent of the excitation energy and geometry of the 
experiment. Secondly, the matrix element Mij includes 
all information concerning the photoemission interaction. 
As such it is sensitive to the experimental geometry via 
symmetry selection rules as well as to the photon-energy. 

The first ARPES experiments on Sr 2 Cu0 2 Cl2 (Ref. 
|l9| ) were able to make the important observation that 
the Zhang- Rice singlet state bandwidth along T — ►(7r, tt) 
is of the order of 2J, and that the Zhang-Rice singlet 
state approaches closest to the chemical potential at the 
(tt/2, 7r/2) point. These observations fit the predictions of 
the t-J-model. However, these early experiments found 
very little dispersion along T — >(tt, 0), which is in con- 
flict with the predictions of the same model. Subsequent 
studies have confirmed the behavior along T— >(7r,7r), 

-Ktt.O)- 



but have also observed dispersion alo ng t he. 
direction in reciprocal space (k-space)c30E2l. 

Nevertheless, the situation as regards the exact disper- 
sion relation along T — >(tt, 0), and in particular concern- 
ing the minimum energy difference between the states 
along r — >(7r, 0) and those at the (tt/2, tt/2) point is un- 
clear. This is more than a mere ARPES detail, as it is 
along the T — >(7r, 0) direction that the predictions of the 
different models vary the most - thus making this direc- 
tion in k-space important for the quantitative comparison 
between theory and experiment. 

Furthermore, the recent controversy surrounding the 
correct Fermi surface topology in the doped hi gh t em- 
perature superconductor Bi 2 Sr2CaCu208+,5 (Ref. 24) has 
illustrated that data sets recorded from the same system 
with different experimental conditions can be remark- 
ably dissimilar. Consequently, both the photon-energy 
and exact polarization geometry used in an ARPES ex- 
periment are important parameters which, if they cannot 
be treated at a quantitative, microscopic level, should at 
least be thoroughly investigated on the experimental side. 

In this paper, we address the electronic structure 
and dynamics of the lowest binding-energy electron- 
removal states in the "standard" undoped model cuprate 
S^CuO^Ck) using ARPES. Following the description of 
the experimental status given above, our "re- visit" of this 
system concentrates on the following points: 

a) a thorough characterization of the photon-energy 
dependence of the first electron-removal states at 
(tt/2, tt/2) and close to (tt/2, 0) (the two points along 
T ^(tt, tt) and T 0) for which the first electron- 
removal states have maximum spectral weight); b) a 
thorough characterization of the polarization dependence 



of the first electron-removal states along these high- 
symmetry directions in k-space; c) the determination of 
the dispersion relation of the first electron-removal states 
along r — >(tt, tt) and T — >(tt, 0), with subsequent compar- 
ison of the results with both existing ARPES data and 
the predictions of an extended t- J-model and d) the de- 
termination of the k-dependent evolution of the coherent 
and incoherent parts of the spectral weight of the first 
electron-removal state along T— >(tt, tt) and r— »(tt, 0), 
and comparison of these data with theoretical predic- 
tions. 



II. EXPERIMENT 

Experiments were performed at the beamlines F2.2 
and W3.2 at the Hamburg Synchrotron Laboratory (HA- 
SYLAB), at the undulator beamline U2-FSGM and 
the 2m SEYA beamline at the Berliner Elcktronen- 
speicherring Gesellschaft fur Synchrotronstrahlung mbH 
(BESSY) . At the storage rings and monochromators used 
in these studies, highly linearly polarized synchrotron ra- 
diation was available. In addition, the crossed undula- 
tor U2-FSGM beamline gave the possibility to use ver- 
tically oriented linearly polarized light. At this latter 
facility, polarization-dependent measurements could be 
performed without changing any other parameter in the 
experiment. A total energy resolution (beamline and an- 
alyzer) of better than 70 meV and an angular acceptance 
of the analyzer of ±1 degree was used. Fresh samples of 
Sr2Cu02Cl2 were cleaved in situ at a base pressure of 
1 x 10 -10 mbar, and spectra were taken within six hours 
after cleavage. Samples were either pre-oriented using 
X-ray diffraction measurements or aligned in-situ with 
the aid of low energy electron diffraction (LEED). In all 
cases, the fine angular adjustment was carried out using 
the k-space symmetry of the sharp ARPES features re- 
lated to non-bonding 02p states around (tt, tt) and (tt, 0), 
as discussed in Refs. |25| and p6| . 

The high-quality single crystalline S^CuO^C^ was 
grown according to the method described in Ref. ||, 
where pre-dried high-purity SrC03, SrCl2 and CuO in 
a ratio 1:1:1 were melted at 1100 °C. Although T N for 
Sr 2 Cu0 2 Cl 2 is 251 K (Ref. |), all ARPES experiments 
were carried out at room temperature. This decision was 
based on two arguments. Firstly, and most importantly, 
as these experiments involve electron ejection from per- 
fect single crystals of a compound with an energy gap 
of the order of 2 eV, we had to eliminate uncertainties 
due to charging effects, which meant measuring at room 
temperature. Secondly, although we measure at a tem- 
perature some 50 K above the Neel temperature, it has 
been shown that the antiferromagnetic spin correlation 
length of Sr 2 Cu0 2 Cl 2 at 350 K is still 250 A (Ref. 
meaning that the lowest binding-energy hole state cre- 
ated in the photoemission process is still embedded in an 
antiferromagnetic spin background. 
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To be absolutely sure that variations in the flux of 
the incident photons from the storage ring did not lead 
to charging-induced energy shifts, we also adjusted the 
beamline such that the photon flux impinging on the sam- 
ple (monitored with a gold mesh upstream of the sam- 
ple) was constant. This means that measuring time is 
the only parameter required for the normalization of the 
data. Bearing in mind that the data presented here are 
a selection of data from about 80 different cleavages, we 
made the qualitative observation that best and sharpest 
energy distribution curves (EDC) were observed for sam- 
ples which behaved most sensitively to charging effects. 
This relation, of course, is reasonable as it correlates 
cleavage quality with the intrinsically highly insulating 
nature of perfect S^CuC^C^. 

Essentially two different experimental geometries were 
used. They can be characterized by the orientation of two 
planes: a) the plane of polarization and b) the emission 
plane. The plane of polarization is the plane spanned by 
the vector of the direction of the synchrotron light and 
the vector of its electric field (vector of polarization). The 
emission plane includes the vector of the direction of the 
photoelectron and the surface normal. Throughout this 
paper we call the experimental geometry parallel if these 
two planes are parallel, perpendicular if these two planes 
are perpendicular to each other (see Fig. || for a sketch 
of the two geometries). 

III. PHOTON-ENERGY AND POLARIZATION 
DEPENDENCE OF THE FIRST 
ELECTRON-REMOVAL STATES OF AN 
UNDOPED CUOa-PLANE 

As mentioned above, both the inconsistencies in the 
literature regarding the dispersion relation of the first 
electron-removal states along the T— >(7r, 0) direction 
(Ref. H|) in Sr 2 Cu02Cl2, as well as the ongoing con- 
troversy regarding the photon-energy dependence of the 
ARPES data from E^S^CaC^Os+s, mean that before 
comparison with theory is carried out, the photon-energy 
dependence of the first electron-removal states should be 
examined in detail. 

In Fig. ||, photon-energy dependent ARPES data for 
the first electron-removal state in Sr2Cu02Cl2 at two 
points in the Brillouin zone: a) (tt/2,tt/2) (perpendic- 
ular geometry); and b) (0.77T, 0) (parallel geometry) are 
shown. In these regions of k-space the spectral weight of 
the first electron-removal state is known to be at-a. lQca l 
maximum along the respective k-space directionst^rEJS. 
The spectra, which are normalized as described in the 
last section, illustrate clearly the strong variation of the 
first electron-removal state intensity with photon-energy. 
The data shown in Fig. ^ represent only a small por- 
tion of the photon-energy dependent data recorded, and 
are intended to give the reader a direct impression of 
the strength of the effects at play. Each of the shown 



spectra is part of a short k|| -series of three or five spec- 
tra. This was done to ensure we always captured the 
spectrum with the highest spectral intensity for the first 
electron-removal states. Within the errors given by the fi- 
nite k-resolution of the experimental setup (0.054 A -1 for 
16 eV photon-energy up to 0.152 A -1 for 80 eV photon- 
energy) we found the highest intensities always at the 
same k-positions, namely at (ir/2,ir/2) and at (0.77T, 0). 

This is in contrast to the results reported in Ref. |23| , 
where series of EDCs on Sr2Cu02Cl2 along T — >(tt^tt) 
taken with two different photon-energies show differences 
not only in intensity, but also in the EDC-derived disper- 
sion relation of the first ionization states. This was dis- 
cussed in Ref. ^3] as not being due to either a) the exper- 
imental setup; b) the excitation of different initial states 
or c) in the variation of kj_ of the photoelectron. The 
effect was rather attributed to the strong impact of the 
matrix element, not only on the strength of the photoc- 
mission signal but also on the position of the quasiparticle 
peak in the photoemission spectra, leading to differences 
in the EDC-derived dispersion relations. We are forced 
to disagree with the last point, as in our extensive col- 
lection of ARPES data, there was never evidence for a 
photon-energy dependent shift of the k-position for which 
the first electron-removal states have minimum binding 
energy along T— >(7r,7r). The same, in fact, holds for 
other directions in the Brillouin zone as well as for the 
dispersion of other features with low binding energy. If 
matrix element effects change the EDC-derived disper- 
sion relation, then this change is, at least for our data, 
smaller than the energy resolution and angular resolving 
power of our ARPES experiments, which is equal to A^ 
or better than A^ the values used in Ref. p3| ). 

Figure [| shows an analysis of the spectral weight of 
the first electron-removal states at the same k-points as 
shown in Fig. ^. These data are derived from a large 
set of photon-energy dependent ARPES data covering 
measurements from 10 (6) cleaves for the upper (lower) 
panels. From Fig. |i| it is clear that the ARPES spectral 
weight of the first electron-removal states oscillates with 
the final state kinetic energy, that is with kj^ . We used a 
value of 8.0 eV for the inner potential E (E = Vq — $, 
where $ is the work function) to calculate k^, a value 
which in the range of 6.9 — 8.9 eV used by other groupscj. 
We see clear maxima at kj_=0.82, 1.63, 2.40 and 3.12 
A -1 , corresponding to photon-energies of 16, 25, 35 and 
48 eV. 

The oscillatory nature of the photon-energy depen- 
dence, coupled to the absence of a classical resonance be- 
havior at the Cu 3p threshold (around 74-76 eV photon- 
energy) indicates that the factor dominating the observed 
behavior is something other than the atomic photoioniza- 
tion cross-sections, and could be related, for example, to 
the extreme two-dimensionality of the electronic states 
concerned. This could lead to a matching of the final 
state kj_ to the periodicity of the unit cell of Si^CuC^C^ 
in c-direction. Note that the differences between the four 
kj^ is 0.72 — 0.81 A -1 which represents in real space a 
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distance of 7.8 — 8.7 A which is comparable to c-axis sep- 
aration of two neighboring CuC>2-planes (= 7.805 A) in 
Sr2Cu02Cl2- The oscillation of the photoemission inten- 
sity with photon-energy can therefore be attributed to 
interference effects of the photoelectron wave diffracted 
from the c-axis periodicity of the layered crystal struc- 
ture, similar to the explanation of the strong photon- 
energy dependence of the photoemission intensity from 
the molecular orbitals of Ceo (Ref. [30|) . 

These strong variations in intensity as a function of 
photon-energy are also reminiscent of ARPES data of 
YBa2Cu307_5, in which an out-of-phase behavior re- 
garding intensity vs. kj_ was observed for the CuC>2- 
dcrived band which crosses the Fermi-energy near the X 
point and the so-called "1 eV peak"Ell. This fact was, at 
that time, used to argue against the surface-state origin 
of the 1 eV peak, a feature which is now believed to rbe 
due to non-mixing 02p states of particular symmetryEi 
Nevertheless, the data presented here, taken together 
with the theoretical and experimental investigations on 
Bi 2 Sr 2 CaCu208+,5 (Refs. [32| and |24|) , indicate clearly that 
care should be taken in the interpretation of absolute 
spectral weights observed in the ARPES data of the lay- 
ered cuprates, as matrix element and diffraction effects 
do play an important role in these systems. 

A further experimental variable in an ARPES ex- 
periment is the polarization of the incoming radiation. 
Figure || shows two series of ARPES measurements on 
S^CuC^Ck) recorded along the high-symmetry direc- 
tions r — >(7r, 0) and r — >(7r, 7r). The photon-energy was 
set to 22.4 eV, which is near a maximum of intensity as 
shown in Fig. || above. In each case, the series are pre- 
sented in pairs of data sets recorded with perpendicular 
and parallel polarization geometries as described in the 
experimental section. 

Along the T— >(tt,tt) direction in k-space the first 
electron-removal state feature shows highest photoemis- 
sion intensity with perpendicular geometry (Fig. |5|a), 
whereas along T — >(tt, 0) maximal intensities are observed 
in the parallel geometry (Fig. ||d), a result which tallies 
with earlier measurements with 25 eV photon-energyo. 
This remarkable dependence of the photoemission inten- 
sity on the polarization was first explained in the context 
of the strong polarization, dependence in photoemission 
data from surface stateso. 

The physical picture behind the polarization depen- 
dence can be described as follows. The interaction op- 
erator A • p has even (+) parity in a parallel and odd 
(-) parity in a perpendicular experimental geometry. As- 
suming the applicability of the Zhang-Rice singlet state 
construction to the first electron-removal final state, in 
a many-body picture this state belongs to the A± g rep- 
resentation, therefore being totally symmetric. Thus in 
this representation the Zhang-Rice singlet has even par- 
ity with respect to a mirror plane along the Cu3d x 2_ y 2- 
02p x /02p y orbital bonds (M 2 ) as well as at 45 degrees 
to the bonds (Mi). In k-space Mi corresponds to the 
r^(7r,7r) direction and to r— >(7r, 0). The initial 



state (ground state) is a one- hole state with d x 2_ y 2- 
symmetry, and therefore has even parity with respect 
to Mi and odd parity with respect to M\. The matrix 
element thus formally vanishes for the two cases 

parallel for T -^(tt, tt) =*> M ~ (+| + |-) = 
perpendicular for T -y(w, 0) M ~ (+| - |+) = 

This argumentation also holds in the one band picture. 
Here the initial state (ground state) is totally symmetric 
and the final state has d x 2_ y 2 -symmetry (as shown in 
Fig. Q), leading to the same result. 

Thus, the observed polarization dependence of the first 
electron-removal states of an undoped Cu02-plane indi- 
cates that these states have a symmetry fully consistent 
with that of the Zhang-Rice singlet. These results amend 
earlier reports (Ref. |2^) regarding the polarization depen- 
dence of the first electron-removal states along V — >(tt, 0) 
in Sr2Cu02Cl2. 

IV. THE DISPERSION RELATION AND 
SPECTRAL FUNCTION OF THE FIRST 
ELECTRON-REMOVAL STATES 

The conupoa picture given by the J-model and its 
extensionsEirtll is a strong dispersion of the Zhang-Rice 
singlet state along the T — >(tt, tt) direction with the min- 
imum binding energy at (7r/2,7r/2). At this k-point, 
the spectral weight of the first electron-removal state 
also has its maximum and vanishes going away from 
(71-/2,71-/2). aepjaLARPES experiments on Sr 2 Cu0 2 Cl 2 
confirm thisEjElaEj. Along the T —>(tt, 0) direction the 
t- J-model predicts a rather low binding energy near the 
high-symmetry point (7r, 0) which is, in fact, almost ener- 
getically degenerate with that at (7r/2, 7r/2). In the same 
model, the spectral weight of the quasiparticle increases 
and is maximal at (tt, 0). This is in contrast to the re- 
sult given in Re f. |l9| where no dispersion was observed 
and in Refs. p0| , [2l[|23| where the binding energy of the 
first electron-removal state increases and the quasiparti- 
cle weight decreases after (7r/2,0). For the (tt, 0)— >(0, tt) 
direction, the t- J-model predicts only little dispersion, 
whereas experiment has shown a strong isotropic dis- 
persion around the (7r/2,7r/2) point (Ref. [H]). Exten- 
sions of the t- J-model and the spin and charge separa- 
tion ansatz exhibit new properties along V — >(7r, 0) and 
(tt, 0)^(0, tt). The Zhang-Rice singlet state around (tt, 0) 
can now lower its energy by derealization and is thus 
pushed to higher binding energies and the dispersion in 
this direction is essentially parabolic. In these models, 
both a quasiparticle dispersion along (tt, 0)— >(0, tt) which 
is isotropic around (tt/2,tt/2) as well as a reduced spec- 
tral weight near (tt, 0) can be achieved. In the spin and 
charge separation model the isotropy along r— >(tt, tt) 
and (tt, 0)— >(0, tt) is intrinsic and near (tt/2, tt/2) the dis- 
persion relation is wedge-like rather than parabolic. Sur- 
prisingly, the data in Ref. [l9| agree well with the extended 
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i- J-model along (tt, 0)— >(0, it), but not along T -^(ir,0). 
In contradiction to that, other data along T — >(tt, 0) are 
consistent with-the extended i- J-model and the spin and 
charge ansatzEu, a v iew w hich appears to be supported 
by further data setsEH'E2l. Thus, the experimental dis- 
persion along the T ^(tt, 0) remains controversial, and 
yet it is the dispersion and the evolution of the spectral 
weight along T —>(tt, 0) that is of the deepest theoretical 
interest for the following reasons. Firstly, the states near 
(-7T, 0) evolve with increasing hole concentration to be- 
come the flat bands located near the chemical potential 
in the metallic systems. These flat band regions are be- 
lieved by some to hold the key to high temperature super- 
conductivity. Secondly, the dispersion and quasiparticle 
spectral weight along T — >(ir, 0) are fingerprints for the 
different models and parameters used in them. Thirdly, 
there are strong indications that the dispersion in the 
insulator is closely connected with the pseudogap in the 
underdoped region of the high T c superconductor phase 
diagram and thus may also be related to the supercon- 
ducting gap in the doped superconducting systemaiU'EE 

Consequently, the importance of experimental data 
which allow the determination of the dispersion and 
spectral weight of the first electron-removal states along 
r — >(tt, 0) cannot be overestimated. Nevertheless, the ex- 
perimentalist is faced with a number of challenges when 
collecting data along this direction: firstly, at several 
points in the Brillouin zone, there is no or only weak 
photointensity from the first electron-removal state. For 
example, near T the Zhang-Rice singlet has no spectral 
weight, because the Gn2d x 2_ y 2 and 02p x /02p y orbitals 
cannot hybridize thereEj, in addition the matrix elements 
for emission along the surface normal are formally zero 
for a perfectly two-dimensional electronic state located 
in the Cu02-plane. Furthermore, near (n, 0) the photoc- 
mission intensity of the first electron-removal states be- 
comes weak as predicted by the extended t- J-models and 
the spin and charge separation ansatz. Secondly, only 
perfectly aligned samples with fresh and clean surfaces 
excited with UV-light in a well-adjusted measurement 
geometry will show usable photocmission data from the 
first electron-removal state. Lastly, as the clean, defect- 
free surface of S^CuO^C^ is highly insulating, avoiding 
charging effects can be difficult. 

In Fig. H we present series of ARPES spectra 
recorded along T —^(tt, 0) taken with parallel polarization 
and r ^(7r,7r) with perpendicular polarization and with 
photon-energies of 16, 22 and 35 eV - thus measuring near 
a peak in the photocmission intensity oscillation in each 
case. The position of the low energy part of the disper- 
sive structure has been fitted using a Gaussian function 
while the high binding-energy background was modelled 
simply with the tail of another Gaussian. Although this 
fit is physically somewhat simplistic, it proved to be suf- 
ficient to reliably find the positions of the peak as indi- 
cated by the triangles in each plot. The absolute value 
on a binding-energy scale changes from one cleavage to 
the next and is rather arbitrary due to the insulating na- 



ture of the substance. However, the two series with 16 
eV were taken during one and the same cleavage, which 
is the prerequisite to compare the binding-energy scales 
of both series. The binding-energy scales of the 22 and 
35 eV data have been adjusted to that of the 16 eV data. 
We present a summary of the dispersion relations derived 
from the fits to the data in Fig. |7[ 

We repeat here that we observe no evidence for a 
photon-energy dependence of the dispersion relation of 
the first ionization states, as is consistent with the ar- 
guments regarding the spectral function given in the in- 
troduction. The dispersion of the first electron-removal 
states along T — >(7r,7r) agrees with the previous results 
having a parabolic shape with its minimum binding en- 
ergy at (n/2, 7r/2). Along T ^(ir, 0) we also find a 
parabolic dispersion, in ag reement with the dispersion 
reported in Refs. ^,|l],^3[ but in contrast to the data 
given in Ref. 19. For 16 eV photon-energy, we measured 
72 meV difference in the lowest binding energies of the 
first electron-removal state along each of the two high- 
symmetry directions, which is half as small as in Ref. |2y . 
In Ref. ^l], a series of spectra representing a "maximum 
intensity cut" for a photon-energy 22.4 eV along a line 
from (ir/2, w/2) to (0.677T, 0) was shown, which displayed 
a dispersion of about 300 meV. However, as we will show, 
the points of maximum intensity along T — >(7r, 0) do not 
coincide with those possessing minimum binding energy. 
Thus such a maximum cut will not trace the line of min- 
imum binding energy. It turns out that this difference 
strongly influences the parameters in model Hamiltoni- 
ans fitted to the dispersion curves. 

We now go on to compare the experimental disper- 
sion with theoretical results obtained within the extended 
t- J-model 



h = -E E + J E s *- s r 



(3) 



1=1 



with the additional hopping terms to second (£2) and 
third (£3) neighbors now added besides the dominating 
nearest neighbor hopping t = t%. The Hamiltonian is 
written in terms of Hubbard operators X[° = c[ a {l — 
ni^ a ) where tr = ±1 is the spin index. The values of 
the hopping terms are determined by mapping the more 
realistic Emery model with parameters derived from a 
constrained density-functional calculationEi to its low 
binding-energy part by means of the cell-perturbation 
method. This procedure gives ^Ai = —0.08 and tz/t\ = 
—0.15. The theoretical dispersion relation shown in Fig. 
is calculated using a variational method involving a spin 
polaron of small radiuses. There-are also related works 
which lead to similar dispersionsEU. 

It has been_shown in exact diagonalisation studies of 
small clustersc3 that the dispersive bandwidth scales with 
the single parameter, J, whereas the form of the disper- 
sion curve is fixed by the ratios tilt\ and t^jt-y. Therefore 
we use the variational ansatz for J = t\ and then we scale 
the bandwidth in Fig. p] with a factor J = 0.22 eV. For 
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reference, we also show the dispersion relation without 
additional hopping terms (ti = ts = 0, scaling factor 
J = 0.28 eV), which shows a much too small energy dif- 
ference between the lowest binding energies of the first 
electron-removal states along V — *(tt, 0) and T —>(tt, it) to 
obtain good agreement with the experimental spectra. 

Another important prediction from model Hamiltoni- 
ans is the evolution of the spectral weight along certain 
directions in the Brillouin zone. Here we map this evo- 
lution and distinguish between the coherent and inco- 
herent parts of the first electron-removal states following 
the procedure proposed in Ref. ^3|, in which the contribu- 
tion from the main valence band tail was first subtracted 
from the spectra. The photoemission intensity was then 
defined as the coherent spectral weight and fitted using a 
Gaussian. The ARPES intensity which is neither part of 
the valence band tail nor in the Gaussian is then taken to 
be a measure of the incoherent part of the first electron- 
removal state spectral function [see Ref. Q . In our case, 
we have used a Lorentzian broadened by the experimental 
resolution (i.e. a Voigt profile) to fit the lowest binding 
energy intensity (the coherent part). Of course, this pro- 
cedure is not physically rigorous, but does offer a rough 
estimation of the possible split between the coherent and 
incoherent spectral weight. 

In Fig. [Sj, the momentum distributions of the two 
parts of the ARPES spectra are shown for 16 eV and 
22 eV excitation energy. Note that the overall shape of 
these distributions is independent of the photon-energy. 
Along r — »(7r, 7r) both the coherent and incoherent parts 
of the first electron-removal states are symmetric around 
(7r/2,7r/2) at which point they both reach their maxi- 
mum weight. In the T — >(ti, 0) direction, we observe a 
steady increase in spectral weight of both components 
up to (0.77T, 0) after which point it drops fast. One can 
find signs of a qualitatively .similar behavior in ARPES 
data from other groupsEjLdSH and for measurements 
on the related substance Ca2Cu02Cl2 (Ref. ^2j), although 
the coherent and incoherent spectral weights were not 
analysed in these studies. Note that along T —*(%, 0) the 
minimum binding energy of the low binding-energy fea- 
ture and the maximum of the coherent intensity do not 
coincide at the same k-point. In calculations usins-ex^ 
tensions of the t- J-model this behavior was predictedE-lo 
and can now be considered as experimentally verified. 

Within the simple fit procedure we can give these quan- 
tities some numbers: the ratio of the coherent to incoher- 
ent spectral weight for the two locations in k-space is 0.6 
((7r/2,7r/2)) and 0.5 ((0.7tt,0)) for 16 eV photon-energy 
and 2.4 ((tt/2, tt/2)) and 3.8 ((0.7tt, 0)) for 22 eV photon- 
energy. It is important to realise that we observe here an 
extremely strong apparent dependence of the ratio be- 
tween the coherent and incoherent spectral weight upon 
the photon-energy. To discuss the physical significance 
of this, we re-write equation ^ as 



S(KE)=- 1 -lm(i\c k ( E _ E 1 (M) _ E 



^inc) c k I *)> 
(4) 



in order to separate the coherent from the incoherent 
spectral weight (£ is the self-energy of the quasiparticlc) . 
Plugging this into equation |l|, it is immediately clear that 
in taking the ratio of the coherent and the incoherent 
intensity, the matrix element cancels out. Thus, this ratio 
is then determined by the spectral function S(k, E) alone 
and should not depend on the excitation energy. 

Our observation that the ratio does appear to depend 
on \w leads then to the following possible conclusions: (i) 
the fit procedure used is too simple and therefore does not 
correctly quantify, however roughly, the coherent and in- 
coherent parts of the spectral function, (ii) there are sig- 
nificant extrinsic contributions to the photoemission sig- 
nal in this binding energy region. These could be the re- 
sult of ao-energetic shift of spectral weight due to inelas- 
tic losseafJ, which could be sensitively dependent on the 
photoelectron kinetic energy, as are the data from elec- 
tron energy loss spectroscopy in reflection of Sr2Cu02Cl2 
(Ref. ^). (iii) the photoemission intensity between the 
the main valence band edge and the chemical potential 
is not derived from the spectral weight (coherent + inco- 
herent) resulting from a single electronic state - i.e. there 
is spectral weight from an additional, different electronic 
state in this energy region. 

The fit procedure, as we have discussed above, is sim- 
plistic - however the coherent-incoherent intensity ratios 
vary by a factor of more than four between the two 
photon-energies. The third possible conclusion - that 
there has to be spectral weight from more than one 
state in this energy region could have important impli- 
cations. We can rule out a significant contribution from 
secondary electrons for the relatively high kinetic ener- 
gies dealt with here—Intensity from surface states is un- 
likely, as XPD datallH shows that the cleavage surface of 
Sr2Cu02Ci2 is terminated within the SrCl units, which 
are essentially ionic and therefore do not support elec- 
tronic states close to the chemical potential. In addi- 
tion, our LEED investigations also gave no evidence for 
a reconstruction of this termination layer. Spin-resolved 
resonant photoemission of CuO (Ref. [!(]) has shown in- 
tensity due to triplet states within leV of the Zhang-Rice 
singlet, although those states would be expected to have 
the same photon-energy dependence as the singlet in our 
experiment. 

If there was an additional, different state in this en- 
ergy region, this would lead to a complicated k- and In- 
dependent overlap between the intensity of this state and 
the higher binding energy components attributed up till 
now to incoherent weight from the Zhang-Rice singlet. 
While this would not be expected to have a large impact 
on the observed dispersion relation for the Zhang-Rice 
singlet (which is, after all, a quantity derived from the 
spectral structure at lowest binding energies for each k- 
point), it could, however lead to the observed photon- 
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energy dependence of the ratio between the low binding- 
energy ("coherent") and higher binding-energy ("inco- 
herent" ) parts of the photoemission spectra. 

V. SUMMARY 

In conclusion, we have presented a detailed ARPES 
study of the low binding-energy occupied electronic 
structure of S^CuC^C^, which corresponds to an inves- 
tigation of the first electron- removal states of an undoped 
Cu02-plane. Our experiments, and the comparison of 
their results with theoretical models, have revealed the 
following main points: 

1. The photoemission signal of the first electron- 
removal states at both (7r/2,7r/2) and (0.77T, 0) exhibits 
a marked photon-energy dependence. The intensity pro- 
file shows strong oscillations with maxima near 16, 25, 
35 and 49 eV, corresponding to final state crystal mo- 
menta k ± =0.82, 1.63, 2.40 and 3.12 A" 1 . This strong 
photon-energy dependence has complicated comparisons 
between data in the literature from different groups as 
regards both the spectral weight and spectral form of 
the first electron-removal states in these systems. We at- 
tribute the oscillation of photoemission intensity (which 
has a period in of ca. 0.75 A -1 ) to the diffraction of 
the photoelectron wave on the periodic c-axis separation 
of the Cu0 2 -planes of 8.2 A. 

2. Along the high-symmetry directions T — >(tt,tt) and 
r — >(7r, 0) the first electron- removal states shows a strong 
polarization dependence. This can be linked to the 
strongly polarization-dependent matrix element, which 
in turn allows the determination of the symmetry of 
the first electron-removal state itself. For both high- 
symmetry directions we observe a polarization depen- 
dence in keeping with that expected for a Zhang-Rice 
singlet state in the framework of either a three-band or 
one-band model Hamiltonian. 

3. Our data show that the dispersion of the first 
electron-removal states along both high symmetry di- 
rections (r — >(7r, 7r) and r — >(7r, 0)) is parabolic-like and 
independent of the excitation energy. This, and the 
rather large difference in lowest binding energy of the 
first electron-removal state along these directions, shows 
the validity of the extended t- J-modcl for describing the 
disperion relation of a single hole in an antiferromag- 
netic Cu02 plane. Thus, the inclusion of second (£2) and 
third (£3) neighbor hopping terms with realistic values of 
t 2 = —0.08 and i 3 = 0.15 in units of the next neighbor 
hopping t = t\ are required. 

4. Upon application of a simple fit procedure, we infer 
the momentum distribution of the spectral weight of the 
coherent and incoherent part of the first electron-removal 
state to have its maximum along T — >(7T, ir) at (ir/2, n/2) 
being symmetrically suppressed away from this point. 
Along r — >(7r, 0) the spectral weights of both parts reach 
their maximum at (0.77T, 0) and then drop fast. The ratio 



between the coherent and incoherent spectral weight is 
strongly photon-energy dependent, which, at first sight 
would appear to violate the physics of the spectral func- 
tion. There are different possible explanations for this in- 
cluding: (i) the necessity for a more sophisticated frame- 
work in which to analyse the weight of the coherent and 
incoherent contributions to the spectral weight; (ii) sig- 
nificant (hi/-dependcnt) intensity due to extrinsic pro- 
cesses (iii) intensity in this energy region due to intrinsic 
electronic states other than the Zhang- Rice singlet. 
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and by the DFG under Fi439/7-l and as part of the 
Graduiertenkolleg "Struktur und Korrelationseffektc im 
Festkorper" der TU Dresden. 
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FIG. 1. Schematic representation of the Zhang-Rice sin- 
glet state in the one-particle representation. The Zhang-Rice 
singlet has the same symmetry as a d a! 2_ !/ 2-orbital, imply- 
ing the existence of two mirror planes perpendicular to the 
Cu02-plane, labelled Mi and M 2 . 



FIG. 2. The experimental geometry for the photoemis- 
sion experiments: a) emission-plane and polarization-plane 
are parallel to each other, b) emission-plane and polariza- 
tion-plane are perpendicular to each other. Throughout this 
paper we call the first case a) parallel and the second case b) 
perpendicular polarization. 
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FIG. 3. A typical example of ARPES data of the first 
electron-removal states of Sr2Cu02Cl2 as a function of pho- 
ton-energy. In panel a) the ARPES EDCs are recorded at 
(7r/2, 7r/2) and in panel b) at (0.77T, 0). The photon-energies 
are as shown. 

FIG. 4. The photon-energy dependence of the photoe- 
mission intensity of the first electron-removal states a) at 
(7r/2, 7r/2) and b) at (0.77r,0), arrived at from the analysis 
of a large body of ARPES data such as that shown in Fig. ^| 
The solid black line represents a guide to the eye and the error 
bars (located in the upper left corner of each panel) indicate 
the error of the fit to the experimental intensity (95% con- 
fidence interval). The top axis shows the corresponding kx 
scale calculated using an inner potential of 8.0 eV (for details 
see text). 



FIG. 8. The coherent and incoherent part of the spectral 
weight of the first electron-removal states determined from the 
data shown in Fig. [jfor a) 22 eV and b) 16 eV photon-energy 
following the procedure used by Leung et al, Phys. Rev. 
B 56, 6320 (1997). The upper part in each case shows a 
representative result of the fit. In the lower part of the figure, 
the weight of the coherent (filled circles) and the incoherent 
(open circles) intensity, normalized to the coherent spectral 
weight at (7r/2,7r/2), is shown. The vertical bars indicate the 
error of the least squares fit to the ARPES intensity (95% 
confidence interval). 



FIG. 5. Series of energy-distribution curves for differ- 
ent polarization geometries. The photon-energy was set to 
22.4 eV and all spectra were taken at 300K. Panels a) and 
b) show measurements along V —>(n,n) where the first elec- 
tron-removal state peak appears near 21.0 eV kinetic energy 
for perpendicular polarization in a). The corresponding se- 
ries along r —>(n,0) are shown in panels c) and d). There, 
the first electron-removal states gives significant photocurrent 
near 20.9 eV kinetic energy only for parallel polarization in 
d). 



FIG. 6. ARPES series showing the dispersion of the 
Zhang-Rice singlet state along the high-symmetry directions. 
The upper three panels show series along V —*(tt, it) the lower 
three panels along V — >(tt, 0). The excitation energy is the 
same for each column, 16 eV (left), 22 eV (center) and 35 
eV (right). Black triangles are included as a guide indicating 
the kinetic energy of the first ionization state determined by 
a fitting procedure described in the text. 



FIG. 7. The dispersion of a single hole in an antiferromag- 
netic Cu02-plane measured along two directions in k-space. 
The symbols represent data from the ARPES spectra from 
Sr2Cu02Cl2, recorded with the indicated photon-energies. 
The vertical error bars in the right panel indicate the errors of 
the fit (95% confidence interval) in determining the binding 
energy of the first ionization states. The k-resolution is 0.054, 
0.075 and 0.094 A" 1 for 16, 22 and 35 eV photon-energy re- 
spectively, which corresponds to 4.5%[6.8%], 6.3%[9.6%] and 
7.8%[11.9%] of the direction T ->(7r,7r)[r ->(vr,0)]. The lines 
indicate the results of the calculations employing a variational 
ansatz using a spin polaron with small radius (for details see 
text). The dashed line corresponds to a calculation where the 
additional transfer terms t2 and are chosen to be zero (i.e. 
the t- J- model), while for the solid line these terms are set to 
t2/t\ — —0.1 and tz/ti = 0.2 (i.e. an extended t- J- model). 
The improved agreement to the experimental data provided 
by the extended f-J-model, particularly along F -^(7r,0) is 
apparent. 



10 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0007283v4 



This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0007283v4 



This figure "fig3.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0007283v4 



This figure "fig4.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0007283v4 



This figure "fig5.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0007283v4 



This figure "fig6.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0007283v4 



This figure "fig7.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0007283v4 



This figure "fig8.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0007283v4 



